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ABSTRACT 

We examine the stellar populations of a sample of 98 z ~ 4.5 Lyman alpha emitting galaxies using 
their broadband colors derived from deep photometry at the MMT. These galaxies were selected by 
narrowband excess from the Large Area Lyman Alpha survey. Twenty-two galaxies are detected in 
two or more of our MMT filters (g', r', i' and z'). By comparing broad and narrowband colors of 
these galaxies to synthetic colors from stellar population models, we determine their ages and stellar 
masses. The highest equivalent width objects have an average age of 4 Myr, consistent with ongoing 
star formation. The lowest EW objects show an age of 40 - 200 Myr, consistent with the expectation 
that larger numbers of stars are causing low EWs. We found masses ranging from 2xl0 7 M Q for 
the youngest objects in the sample to 2xl0 9 M Q for the oldest. It is possible that dust effects could 
produce large equivalent widths even in older populations by allowing the Lyman alpha photons to 
escape, even while the continuum is extinguished, and we present models for this scenario also. 
Subject headings: galaxies: fundamental parameters - galaxies: high-redshift - galaxies: ISM - galax- 
ies: evolution 



1. INTRODUCTION 

There are two popular techniques for locating galax- 
ies at high redsh ift: the Lyman Break technique 
(|Steidel et alJ H~996) and narrowband selection. The sec- 
ond entails observing a galaxy in a narrowband filter con- 
taining a rcdshifted emission line. This method is useful 
to select high redshift galaxies with strong Lyman alpha 
(Lya) emission lines. While the line emission proper- 
ties of Lyman Break galaxies (LBGs) are known, little 
is known about the continuum properties of Lya galax- 
ies because they are so faint. Studying the continuum 
properties of these Lya galaxies at high-redshift can tell 
us about their age, stellar mass, dust content and star 
formation history. 

We have used the narrowband selection technique to lo- 
cate and study Lyman alpha emitting galaxies (LAEs) at 
z ~ 4.5. Similar studies have been done both in this field 
and in other fields (e.g., Rhoads et al 2000, 2004; Rhoads 
& Malhotra 2001; Malhotra & Rhoads 2002; Cowie & Hu 
1998; Hu et al 1998, 2002, 2004; Kudritzki et al 2000; 
Fynbo, Moller, & Thomsen 2001; Pentericci et al 2000; 
Stiavelli et al 2001; Ouchi et al 2001, 2003, 2004; Fu- 
jita et al 2003; Shimasaku et al 2003, 2006; Kodaira et 
al 2003; Ajiki et al 2004; Taniguchi et al 2005; Vene- 
mans et al 2002, 2004). In many of these studies, the 
strength of the Lyman alpha (Lya) line has been found 
to be much greater than that of a normal stellar pop- 
ulation (Kudritzki et al. 2000; Dawson et al. 2004). 
Malhotra and Rhoads (2002; hereafter MR02) found nu- 
merous LAEs with rest frame equivalent widths (EWs) > 
200 A. Assuming a constant star formation rate, the EW 
of a normal stellar population will asymptote toward 80 
A by 10 8 years, down from a maximum value of ~ 240 
A. These galaxies are of interest because it is believed 
that this strong Lya emission is a sign of ongoing star 



form ation activity in these galaxies ([Partridge fc Peebles! 
Il967f ). and it could mean that these are some of the 
youngest galaxies in the early universe. 

There are a few possible scenarios that could be cre- 
ating this large EW. Strong Lya emission could be pro- 
duced via star formation if the stellar photospheres are 
hotter than normal, which could happen in low metallic- 
ity galaxies. It could also mean that these galaxies have 
their stellar mass distributed via a top-heavy initial mass 
function (IMF). Both of these scenarios are possible in 
primitive galaxies, which are thought to contain young 
stars and little dust. 

Active galactic nuclei also produce a large Lya EW. In 
an optical spectroscopic survey of z ~ 4.5 Lya galaxies, 
Dawson et al. (2004) found large Lya equivalent widths, 
but narrow physical widths (Av < 500 km s~ 4 ). They 
also placed tight upper limits on the flux of accompany- 
ing high-ionization state emission lines (e.g., NV A1240, 
SilV A1398, CIV A1549 and Hell A1640), suggesting that 
the large Lya EW is powered by star formation rather 
than AGN. Additionally, Type I (broad-lined) AGN are 
ruled out because the width of the Lya line in such an 
object is greater than the width of the narrowband filters 
which were used (see Section 2.3). While Type II AGN 
lines are narrower than our filters, they remain broader 
than the typical line width of our Lya sample as de- 
termined by followup spectroscopy (e.g., Rhoads et al. 
2003). Malhotra et al. (2003) and Wang et al. (2004) 
searched for a correlation between LAEs and Type II 
AGN using deep Chandra X-Ray Observatory images. 
101 Lya emitters were known to lie in the two fields 
observed with Chandra, and none of them were detected 
at the 3er level (L 2 - S keV = 2.8 x 10 42 ergs s" 1 ). The 
sources remained undetected when they were stacked, 
and they concluded that less than 4.8% of the Lya emit- 
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ters they studied could be possible AGNs. Therefore, we 
are confident that there are at most a few AGN in our 
Lya galaxy sample. 

Another possible scenario involves enhancement of the 
Lya line in a clumpy medium. In this scenario, the con- 
tinuum light is attenuated by dust whereas the Lya line 
is not, resulting in a greatly enhanced observed Lya EW 
(Neufeld 1991; Hansen & Oh 2006). This can happen if 
the dust is primarily in cold, neutral clouds, whereas the 
inter-cloud medium is hot and mainly ionized. Because 
Lya photons are resonantly scattered, they are preferen- 
tially absorbed at the surface. Thus it is highly unlikely 
that a Lya photon will penetrate deep into one of these 
clouds. They will be absorbed and re-emitted right at the 
surface, effectively scattering off of the clouds and spend- 
ing the majority of their time in the inter-cloud medium. 
However, continuum photons are not resonantly scat- 
tered, and if the covering factor of clouds is large, they 
will in general pass through the interior of one. Thus the 
continuum photons will suffer greater attenuation than 
the Lya photons, effectively enhancing the Lya EW. 

These Lya galaxies are so faint compared to other 
high-z objects that have been studied, that not a lot is 
known about their continuum properties. Our primary 
goal is to better constrain the ages and stellar masses of 
these objects from their continuum light in order to bet- 
ter understand the driving force behind the strength of 
the Lya line. We have obtained deep broadband images 
in the g' , r', i' and z' bands for this purpose. Using the 
colors of these galaxies, we can study continuum proper- 
ties of individual Lya emitting galaxies at this redshift 
for the first time, allowing us to estimate their age, mass, 
star formation rate and dust content. The next phase of 
our project will be to look into the likelihood of dust 
causing the large Lya EW, which we will begin to do 
in this paper, as the colors of these galaxies might dis- 
tinguish the cause of the large EW. Blue colors would 
indicate young stars with hot photospheres, while red 
colors would indicate dust quenching of the continuum, 
enhancing the Lya EW. Observations, data reduction 
technique and sample selection are reported in Section 
2; results are presented in Section 3, and they are dis- 
cussed in Section 4. Conclusions are presented in Section 
5. 

2. DATA HANDLING 

2.1. Observations 

The Large Area Lyman Alpha (LALA) survey began in 
1998 using the 4m Mayall telescope at the Kitt Peak Na- 
tional Observatory (KPNO). The final area of this survey 
was 0.72 deg 2 in two fields, Bootes and Cetus, centered 
at 14 h 25 m 57 s , +35 Q 3 2' (J2000) and 02^0 5 m 20 s , -04°55' 
(J2000) respectively (jRhoads et all l2000r i. The LALA 
survey found large samples of Lya emitting galaxies at 
z = 4.5, 5.7 and 6.5, with a spectroscopic success rate 
of up to 70%. We observed the LALA Cetus Field for 
three full nights in 2005 November, and again for four 
1/4 ni ghts in 2006 Januar y, using the Megacam instru- 
ment (M cLeod et~allll998| l at the Multiple Mirror Tele- 
scope (MMT). Megacam is a large mosaic CCD camera 
with a 24'x24' field of view, made up of 36 CCDs with 
2048x4608 pixels. Each pixel is Of! 08 on the sky. We ac- 
quired deep broadband images in four Sloan Digital Sky 
Survey (SDSS) filters: g', r', i' and z'. The total expo- 



sure time in the g' , r 1 , i' and z 1 filters were 4.33, 3.50, 
4.78 and 5.33 hours respectively. Given that the seeing 
during the run was rarely below 0'.'8, we chose to bin the 
pixels 2x2 to reduce data volume, resulting in a final 
pixel scale of ~0'/16 per pixel in our individual images. 

2.2. Data Reduction 

The data were red uced in IRAF 1 (lTodylll986L I1993J) 
using the MSCRE D dValdes fc Todvlll998t lValdelfl998h 
and MEGARED (jMcLeod et all l2006h reduction pack- 
ages. The data were first processed using the standard 
CCD reduction steps done by the task cedproc (overscan 
correction, trim, bias subtraction and flat fielding). Bad 
pixels were flagged using the Megacam bad pixel masks 
included with the MEGARED package. 

Fringing was present in the i' and z' band data, and 
needed to be removed. For this, we made object masks 
using the task objmasks 2 . We then took a set of sci- 
ence frames in each filter, and combined them to make a 
skyflat using sflatcombine (using the object masks to ex- 
clude any objects from the resultant image). This skyflat 
was then median smoothed on a scale of 150x150 pixels 
(i') and 200x200 pixels (z 1 ). The smoothed images were 
subtracted from the unsmoothed images, and the result 
was a fringe frame for each band. The fringing was then 
removed from the science images using the task rmfringe, 
which scales and subtracts the fringe frame. To make 
the illumination correction, skyflats were made, median 
smoothed 5x5 pixels, normalized and divided out of the 
science images. 

The amplifiers in the reduced images were merged 
(from 72 to 36 extensions) using the MEGARED task 
megamerge. The WCS written in the image header at 
the telescope is systematically off by 0.5 degrees in rota- 
tion, so we used the task fixmosaic to adjust the rotation 
angle so that the WCS solution works better. The WCS 
was then determined using megawes, and the WCS dis- 
tortion terms were installed in the header using the task 
zpn. 

The 36 chips were combined into a single extension 
using the program SWarp 3 . We ran SWarp on each in- 
put image, using a corresponding weight map which gave 
zero weight to known bad pixels and cosmic ray hits. 
We located the cosmic rays using the method of Rhoads 
(2000). We specified the center of the output images to 
be the average center of the five dither positions we used 
(02 ft 04 s 56 m , -05°01'01" J2000), which corresponded to 
an image size of 11988x10881 pixels and a pixel scale of 
0'.'1587 per pixel. During this process, each image was 
resampled (using the Lanczos-3 6-tap filter) and inter- 
polated onto the new pixel grid. While this was done, 
SWarp determined and subtracted the background value 
in each extension. 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories (NOAO), which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 

2 In order to prevent fringing from making it onto the object 
masks, we first made a skyflat without object masks, and divided 
this skyflat out of the science data. We then got the object masks 
from these quasi-fringe removed images. The rest of the procedure 
used the untouched science frames. 

3 SWarp is a program that resamples and co-adds together FITS 
images, authored by Emmanuel Bertin. It is distributed by Terapix 
at: http: / /terapix.iap.fr/soft/swarp 
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An initial stack for each filter was made using SWarp 
to average together the individual images. To make the 
final version of the stack, we needed to remove the satel- 
lite trails which littered our images. To do this, we cre- 
ated a difference image between each individual input 
image and the initial stack. Using a thresholding tech- 
nique, we flagged the satellite trails in this image, and 
combined those flags with the existing weight-maps. We 
ran SWarp a final time with this as the input weight 
map, creating the final stack. We did this process for 
each band, giving us four images, one final stack in each 
of the g' , r', i 1 and z 1 bands. In an attempt to increase 
the signal of our objects, we tried co-adding R and I 
band imag es from the NOAO D eep Wide-Field Survey 
(NDWFS) (jJannuzi fc Devlll999f) to our r' and i' images, 
however the additional signal did not make a significant 
difference in the results. 

To calculate the photometric zeropoint for our obser- 
vations we took images of the standard stars SA95 190 
and SA95 193 during our run, where the zeropoint is 
defined as the magnitude of an object with one count 
in the image (for an integration time of 400s). These 
standard stars are from Landolt (1992), and we used the 
transforms from Johnson-Morgan-Cousins to SDSS mag- 
nitudes from Fukugita et al. (1996). Averaging the re- 
sults from the two stars for each band, we obtained these 
zeropoints: #'=33. 13, r'=33.00, i'=32.58 and z'=31.43. 
Using these zeropoints with our data, all of our results 
are in AB magnitudes 4 . 

2.3. Lyman a Galaxy Selection 

To extract the ob jects from each stack, w e used the 
SExtractor package (jBertin &: Arnouts1 [l996). For SEx- 
tractor to accurately estimate the errors in the flux, it 
needed to know the gain in the input images. To calcu- 
late this, we needed the mean and standard deviation of 
the background of the images. Because the background 
was subtracted out in SWarp, we obtained the average 
background value in each band by averaging the values 
obtained from the pre-SWarped images, and adding this 
value onto the final stacks. We then ran SExtractor in the 
two-image mode using a nine pixel aperture (2. "32) with 
the narrowband images from the LALA survey, giving us 
catalogs of objects which were detected in both images. 

The method we have used to locate Lya galaxies in- 
volves taking an image using a narrowband filter con- 
taining the wavelength for Lya at a certain redshift. We 
used the catalog from MR02, which was based on LALA 
narrowband images and broad band B Wl R, and / data 
from the NDWFS (Jannuzi & Dey 1999). The narrow- 
band data consists of five overlapping narrowband filters, 
each with a FWHM ~ 80A. The central wavelengths are 
AA6559 (HO), 6611 (H4), 6650 (H8), 6692 (H12) and 6730 
(H16), giving a total redshift coverage 4.37 < z < 4.57. 
To ensure that there was no overlap, we only used objects 
selected from the HO, H8 and H16 filters in our analy- 
sis. In order to compare the two data sets, the Megacam 
data were registered and remapped onto the same scale 
as the narrowband data using the IRAF tasks wesmap 
and geotran. 

Selection criteria for the MR02 catalog, following 



Rhoads & Malhotra 2001, were as follows: (1) 5<r signif- 
icance detection in the narrowband: This is calculated 
using an SExtractor aperture flux with the associated 
flux error; flux/error > 5. (2) 4er significant excess of 
narrowband flux: This was calculated by taking narrow- 
band flux - broadband flux, both calibrated in physi- 
cal units, and the associated error (sqrt(error 2 b roac ; + 
error 2 narrcra ,)) and demanding that flux difference / er- 
ror > 4. (3) Factor of > 2 ratio between broad and 
narrow band fluxes, calibrated to the same units. (4) No 
more than 2 a significant flux in the bluest filter observed 
(B w ). 

2.4. Stellar Population Models 

In order to study the properties of the galaxies in 
our sample we compare them to stellar population mod- 
els, using the stellar population modeling software by 
Bruzual & Chariot (2003). With these models, we were 
able to choose a range of ages, metallicities and star for- 
mation rates for comparison with our sample. We chose 
ages ranging from 10 6 - 10 9 years, metallicities from .02 
Z Q - Z Q and exponentially decaying star formation rates 
with a characteristic time-scale of r = 10 3 - 2x 10 9 years. 
We included dust via the Calzetti dust extinction law 
(Calzetti et al. 1994), which is applicable to starburst 
galaxies, in the range: A1200 = 0-2. Lastly, we included 
intergalactic medium (IGM) absorption via the prescrip- 
tion of Madau (1995). 

The Bruzual & Chariot code (BC03) output the flux of 
a given stellar population in units of L Q A -1 . In order to 
go from flux in these units to bandpass averaged fluxes, 
we used the method outlined by Papovich et al. (2001). 
In short, we took the output from BC03, and converted 
it from i\ into i v (units of erg s _1 cm~ 2 Hz -1 ). This 
flux was then multiplied by the transmission function for 
a given bandpass, and integrated over all frequency. The 
bandpass averaged flux (/„) is this result normalized to 
the in tegral of the transm ission function. AB magni- 
tudes ()Oke h Gunnlfl983T l for the candidates were then 
computed. 

In order to calculate colors which we could directly 
compare to our observations, it was necessary to add in 
emission line flux from the Lya line, which appears in 
the r' filter at z ~ 4.5. While the BC03 software does 
not calculate emission line strengths, we were able to 
use one of its many other output quantities: the number 
of ionizing photons. We calculated the Lya emission line 
strength (in units of Lq A -1 to match the model output) 
by using: 

he 2 

Lya Linestrength = - — — x — x n; lon (1) 

ALyaLQtXA 3 

where iii on is the number of ionizing photons, A A is the 
bin size of the wavelength array (1A), and the factor 
of 2/3 represents the fraction of ionizing photons which 
will produce Lya photons when interacting with the local 
interstellar medium (ISM) under Case B recombination 5 . 
In order to model the clumpy dust scenario, we needed to 
ensure that the Lya line did not suffer dust attenuation. 
To do this, the continuum flux was multiplied by the 
Calzetti dust law before we added in the Lya flux to the 



For more details on our red uction and analysis process, see: 
http://stcvcnf.asu.edu/Reduction.html 



5 This calculation follows the simple assumptions that no ioniz- 
ing photons escape, and that all Lya photons escape. 
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spectrum at the correct wavelength bin. In this case, the 
continuum suffers dust attenuation while the Lya line 
does not. 

3. RESULTS 
3.1. Lya Galaxy Candidates 

Using the catalog described above, we have identified 
98 objects as Lya galaxy candidates within the 24'x24' 
Megacam field of view. We have classified 22 of these 
objects as continuum detections on the basis that they 
have at least 2a detections in two of the r', i' or z' bands 
(the other 76 objects were undetected at this level in the 
broadband data, but we stacked their fluxes for analysis). 

We have in our possession IMACS spectra of objects in 
the LALA Cetus Field (Wang et al. 2007). These spec- 
tra have been used to identify the redshift of the Lya line 
in these galaxies. Out of the 98 total Lya galaxy candi- 
dates, 28 have been spectroscopically confirmed as being 
Lya emitters at z ~ 4.5. Seven of these confirmations 
are among the 22 Megacam detections. Another seven 
of the Megacam detections have IMACS spectra but had 
no strong Lya line identified, however some of these may 
show a line with deeper spectra. 

3.2. Equivalent Width Distribution 

We calculated the rest-frame equivalent widths of our 
22 detections using the ratio of the line flux to the con- 
tinuum flux via: 

EW = ( n 1-77 j ) x (1 + z)" 1 (2) 

\ AA BB AAnb / 

where r\ is the ratio of narrowband to broadband flux (in 
units of f„), and the A As are the filter widths (80A and 
1220A for the narrow and broad filters respectively). The 
narrowband fluxes are measured from the LALA narrow- 
band data discussed above, and the broadband fluxes are 
measured from our MMT r' data. 

Figure 1 shows the calculated rest-frame EWs of our 
22 continuum detected galaxies. The range of EWs is 5 — 
800 A 6 . While the range of our EW distribution peaks at 
the low end of previous studies (i.e. most of our objects 
have rest- frame EW < 200A), we recognize that the 76 
undetected objects would all lie at the high EW end of 
this plot (in fact, they peak at EW well over 200 A). Be- 
cause they are undetected, this means that their broad- 
band fluxes were too faint to be detected in our survey. 
However, they still have very bright narrowband fluxes, 
indicating strong Lya emission. This strong emission, 
coupled with faint broadband r' magnitudes, results in a 
high Lya EW. 

From models of normal stellar populations (those with 
a near-Solar mctallicity and a Salpeter IMF), it was 
found that the maximum EW a galaxy could have is 
~ 240A at t ~ 10 6 yr, with the EW approaching a 
steady value of 80A at t - 10 8 yr (MR02; Chariot & 
Fall 1993). Our data show many objects with a rest- 
frame EW greater than 100 A. While individual galaxies 

6 However, the object with an EW of 5 A is likely an interloper. 
It has an ~ 10 <r detection in g' , and > 30 a detections in r', i'and 
z' . We will exclude this object from the rest of our analysis, giving 
us 21 detections. 



with EW this high are possible from normal stellar popu- 
lations, the ratio of high-low EWs is much greater than it 
ought to be when we account for the 76 objects without 
MMT detections. This implies that something is causing 
the EW to be higher than normal. One of our goals in 
this study is to investigate whether the cause is massive 
stars creating more Lya photons, or dust suppressing the 
continuum and thus enhancing the Lya EW. 
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Fig. 1. — The distribution of rest-frame Lya equivalent widths 
from our 22 detections. The EWs range from 5 - 800 A. The single 
object with EW < 20 A is likely an interloper that has managed 
to satisfy our selection criteria. A normal stellar population has an 
EW that asymptotes toward ~ 80 A as its age goes to 10 8 years, 
so something is causing the Lya EW of many of these objects to 
be higher.. 

We note that four objects among our MMT detections 
show observed equivalent width in excess of 230A. Such 
high equivalent widths might be explained by a combina- 
tion of extreme youth and/or low metallicity (e.g. Mal- 
hotra & Rhoads 2004), or might be produced by radia- 
tive transfer effects (Neufeld 1991, Hansen & Oh 2006). 
However, two of these objects are consistent with 200A 
equivalent width at the la level, and the other two at the 
2a level, so such mechanisms are not strictly required to 
explain these objects. Stronger evidence for such effects 
may be found among the 76 sources whose continuum 
emission was too faint for our MMT data. 

3.3. Stacking Analysis 

Our goal has been to be able to compare our objects 
to the models to obtain estimates of physical parameters 
such as age and stellar mass, along with dust content. 
We divided our objects into six groups which we then 
stacked (in order to damp down the errors) to get com- 
posite fluxes which we could compare to the models. The 
groups arc: non-detections, which consists of the 76 se- 
lected objects which did not have a 2 a detection in at 
least two bands (out of r', i' and z'); low EW, which 
consists of the seven detected objects with EW of 20- 
40A; middle EW, which consists of the seven detected 
objects with EW of 45-100A; high EW, which consists of 
the seven detected objects with EW > llOA; detections, 
which consists of the 21 objects which were detected in 
the Megacam data; spectroscopically confirmed, which 
consists of the 28 objects (7 detected) which were con- 
firmed to be Lya emitters at z ~ 4.5. These stacks are 
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shown in Figure 2. 

3.4. Age and Mass Estimates 

In order to study these galaxies, we have created nu- 
merous synthetic stellar population models which we 
compare to our observed galaxies. The most illuminat- 
ing way in which we can study the observed vs. model 
galaxies are in a color-color plot (Figures 2 and 3). We 
have plotted the r' - i' color vs. the r' - nb color of our 
objects, and then overplotted many model curves. Be- 
cause changing the metallicity did not much change the 
position of the model curves in the plane we are study- 
ing, we have elected to only use models with Z = .02 
Z Q . The differing line styles represent the two different 
star formation rates we used, with the solid lines repre- 
senting continuous star formation, and the dashed lines 
representing exponentially decaying star formation with 
a decay time of r = 10 7 years. 

We ran the models at 24 different ages ranging from 
1 Myr - 2 Gyr 7 , and in the figure we show the five ages 
which best surround the data, represented by different 
colors. We included dust in the models, ranging from 
A1200 = 0-2. The extent of the model tracks repre- 
sents the different dust optical depths. As we discussed 
above, this dust does not affect the strength of the model 
Lya line. We have connected the zero-dust end of the 
model tracks with the dotted black line in order to show 
where galaxies with no dust should lie. 
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Fig. 2. — Color-color plot showing the six stacked points along 
with stellar population model tracks. The vertical axis is the r' - 
narrowband (HO) color, and the horizontal axis is the r' - i' color. 
The model tracks shown represent the colors of models in this 
plane. All models are for .02 7,q . The solid tracks denote contin- 
uous star formation, while the dashed tracks denote exponentially 
decaying star formation, with a decay time of r = 10 7 years. The 
different colors of the model curves represent the ages of the mod- 
els. The length of the model curves represent the model colors with 
differing amounts of dust. Models were run with A1200 =0 — 2, 
with the A1200 point lying at the smallest r' - nb color. The dotted 
black line connects the bottom of the model curves, and as such is 
the "zero-dust" line. Ages were found for each stack by minimizing 
the distance from the stack to the nearest A1200 = model point 
(i.e. nearest point on the zero-dust line) for each star formation 
rate (20 different possible ages were used). 



With significant flux in only three broadband points 



(r' , i' and z'), it would be hard to fit each stack to mod- 
els with a full range of every parameter. However, given 
what we have learned from studying these objects in the 
color-color plane, it is relatively straightforward to fit an 
age to each of these stacks (assuming a star formation 
history and no dust). In order to fit the age, we con- 
sidered 20 different ages ranging from 1 Myr - 200 Myr 
(including the five ages shown in Figure 2). By finding 
the model point (we used the A1200 = point for each 
model; see discussion for details) closest to each stack, we 
have effectively found the average age of an object in the 
stack. Because we didn't want to limit ourselves to one 
star formation history, we have found two ages for each 
object; one from the closest model with a constant star 
formation rate (SFR), and one from the closest model 
with an exponential SFR. These ages are reported in Ta- 
bles 1 and 2 for constant and exponentially decaying SFR 
respectively. For a constant SFR, the ages ranged from 
1 - 200 Myr, while for an exponentially decaying SFR, 
they ranged from 1-40 Myr. 

In order to find the stellar mass, we wanted to only 
use the best-fit model for each stack. We took that to 
be the model with the best-fit age found above, and we 
calculated a mass for each SFR. The mass was found by 
a simple weighted ratio of object flux to model flux. We 
derived the weighting by minimizing the ratio of the mass 
error to the mass. The mass error was calculated via a 
monte carlo simulation using the flux errors from each 
band for each stack. These mass estimates are reported 
in Tables 1 and 2. The masses of the objects from both 
types of star formation rates range from-2xl0 7 -2xl0 9 
M Q . These masses are indicative of an average object of 
the stack, and likewise the reported error is the error in 
the mean of the object mass in each stack, rather than 
the error in the measurement. 




Fig. 3. — Color-color plot showing the location of the 21 indi- 
vidual detections (along with the outlier), along with stellar pop- 
ulation model tracks. The individual galaxies are plotted with 
different colored symbols to denote whether they have been spec- 
troscopically confirmed to be at z ~ 4.5. The axes and model tracks 
are the same as in Figure 2. In general, objects above the zero- 
dust line would have their EW enhanced due to dust effects, while 
objects below this line would have intrinsically high EWs, or have 
their whole spectrum attenuated by a geometrically homogeneous 
dust distribution (for the low EW subsample). 



7 At 1 Myr, the two different star formation rates have identical 
colors, so their model tracks overlap, that is why it appears as if 
there is no curve for an exponentially decaying galaxy at 1 Myr. 



In Figure 3, we have plotted the 21 individually de- 
tected galaxies (plus the one outlier), along with their 
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TABLE 1 

Best-Fit Ages and Masses for a Constant SFR 



Stack 


r' mag. 


Age 


x 2 


Mass (1O 8 M ) 




of stack 


(Myr) 




per object 


7 Low EW 


23.27±0.10 


200 


18.48 


23.54±1.66 


7 Mid EW 


23.49±0.12 


80 


0.93 


8.52±0.74 


7 High EW 


24.08±0.21 


4 


0.94 


0.95±0.15 


21 Dct. 


22.39±0.09 


40 


0.62 


4.60±0.29 


76 Undet. 


23.31±0.47 


1 


6.59 


0.19±0.07 


28 Spec. Conf. 


22.96±0.18 


3 


0.05 


0.62±0.09 



Note. — Each model has a metallicity of .02 Z Q . The 
low EW bin is from 20-40 A; the middle bin is from 45-100 
A; the highest bin has objects with EW > 110 A. The bin 
with 21 detected objects consists of the 21 galaxies which 
were detected at the 2a level in at least two of the r', i' or 
z' bands. The last bin consists of the 28 objects which were 
spectroscopically confirmed to have a Lya line at z ~ 4.5. 
The x 2 is defined as the square of the distance in sigma units 
to the best-fit age point (See Figure 2). The r' magnitude is 
for the stack as a whole, while the mass is per object. 

la error bars. We have distinguished between the indi- 
vidual galaxies in order to show those that were spectro- 
scopically confirmed. By their position in this plot, we 
hope to find out whether they have an intrinsically strong 
Lya line, or one that is enhanced due to dust effects. 

4. DISCUSSION 
4.1. Age and Mass Estimates 

We fitted the ages of our stacked points using dust-free 
models. Objects in the stack of the 76 non-detections 
have the lowest age by far of all of the stacks (1 Myr), 
and thus the lowest mass (0.19 x 10 8 Mq ). This is not 
surprising, however, because these objects are selected 
based on their Lya line strength. In these galaxies, the 
line is strong enough to be seen clearly but the continuum 
is not, so they are intrinsically faint galaxies, implying a 
low stellar mass. Their colors are best fit by a low age, 
which also points to a low mass; the galaxy has simply 
not had that much time to form stars. 

The three equivalent width bins show a logical trend 
in age and mass. The lowest EW bin has the highest 
age at 200 (40) Myr, and thus the highest mass at 23.54 
(16.15) x 10 8 Mq for a continuous (exponentially decay- 
ing) SFR. This age is consistent with the fact that the 
EWs in this bin were not higher than normal, and thus 
this EW bin is best fit by an older stellar population. 
The middle EW bin has a younger age and a smaller 
mass, and the highest EW bin has the youngest age out 
of the three bins at 4 (4) Myr, with the smallest mass at 
0.95 (0.94) x 10 8 M Q . This would indicate that the ob- 
jects with the highest EW are truly the most young and 
primitive galaxies in our sample (that were detected). 
The results for the stack of all 21 MMT-detected objects 
are intermediate among the results grouped by equivalent 
width, as expected. The ages for this stack were 40 (10) 
Myr, and masses were 4.6 (2.0) x 10 8 M . The stack of 
spectroscopically confirmed objects should lie somewhere 
between the detected and non-detected stacks, and in- 
deed this is the case, with ages of 3 (2) Myr, and masses 
of 0.62 (0.68) x 10 8 Mq . 

The overall mass range we found of ~ 2xl0 7 - 2xl0 9 





TABLE 2 






Best-Fit Ages and Masses for an Exponentially Decaying 




SFR 






Stack 


Age 


x 2 


Mass (10 8 M Q ) 




(Myr) 




per object 


7 Low EW Objects 


40 


4.05 


16.15±1.14 


7 Mid EW Objects 


20 


1.05 


4.24±0.37 


7 High EW Objects 


4 


0.97 


0.94±0.14 


21 Det. Objects 


10 


3.43 


1.97±0.13 


76 Undet. Objects 


1 


6.59 


0.19±0.07 


28 Spec. Conf. Objects 


2 


0.09 


0.68±0.10 



Note. — See Table 1 notes. 

Mq is comparable to other studies of similar objects. Due 
to the fact that our observed magnitudes are brighter 
than the majority of surveys which we can compare to, 
we point to the trend between apparent magnitude and 
mass. That is, that the fainter an object is, the less 
massive it ought to be (assuming they are at similar red- 
shifts). Ellis et al. (2001) discovered a very faint (I ~ 
30) Lya emitting galaxy at z = 5.576 that was able to be 
detected because it was lensed by a foreground galaxy. 
Due to its faint magnitude, this object was determined 
to have a mass of around 10 6 Mq , implying a very young 
age (~ 2 Myr). This object is consistent with our results, 
in that it is both much fainter and much less massive than 
the objects we are studying in this paper. Gawiser et al. 
(2006) studied numerous LAEs which they detected at 
z^ 3.1, with a median R magnitude of ~ 27. In order 
to take advantage of the full spectral range available to 
them, they stacked their sources before performing SED 
fitting. When this was done, they found an average mass 
per object of ~ 5 x 10 s Mq . All of our stacks (except the 
stack of undetected objects) have an average individual 
object r' magnitude of less than 27, so even given the dif- 
ferences in photometric systems, they are likely brighter 
(much brighter in fact, due to their increased distance). 
While (depending on the SFR) not all of these stacks 
have masses higher than 5 x 10 8 M Q , they are of the 
same order. 

Perhaps the best comparison to our results comes from 
the study by Pirzkal et al.(2007), where they are study- 
ing the properties of z ~ 5 LAEs detected in the GRism 
ACS Program for Extra galactic Science (GRAPES) sur- 
vey (jPirzkal et alJ[200 4) . Using models with an exponen- 
tially decaying SFR, they found a mass range of 3 x 10 6 
- 3 x 10 s M Q in objects with an i' range of 25.52 - 29.35 
(the i' range of an average object in our stacks is from 
25.12 - 29.90). While the majority of our objects have 
i' brighter than 26 (all except objects from the unde- 
tected stack and the spectroscopically confirmed stack), 
all of the GRAPES objects except one have i' fainter 
than 26. This explains well the fact that our derived 
masses are in general greater than those derived for the 
GRAPES LAEs. We are studying intrinsically brighter, 
and thus more massive objects. While our wavelength 
baseline is not very large, we are able to estimate the 
mass of stars producing UV and Lyalight. 

4.2. Dusty Scenario 

In studying Figure 3, we have identified three distinct 
regions in the color-color plane: one region which would 
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house intrinsically blue objects, one region of red objects 
with high EW, and one of red objects with low EW. 
The region that lies below the zero-dust line and blue- 
ward of r' - i' = would appear to house objects that 
have moderate-to-high EWs with blue colors, indicating 
that they are intrinsically blue, containing young stars. 
The lower EW/redder r' - i' color sub- area of this region 
would contain objects may have some dust, but it is af- 
fecting both the continuum and Lya line, lowering the 
EW and reddening the color. 

The second region consists of the area above the zero- 
dust line with an r' - nb color > 2. This region would 
hold objects that lie in the dusty regime of the models, 
yet they would still have "red" r' - nb colors, indicating 
a large EW. Objects such as these could be explained 
by dust quenching of the continuum (causing the red 
r' - i' color), while the Lya line would not be affected 
by the dust, making it appear to be strong (causing the 
red r' - nb color). The last region consists of the area 
above the zero dust line, but with a smaller r' - nb color 
(r' - i' ~ 1, r' - nb = 1 - 1.5). Objects in this area 
would have a redder r' - i' color, with a smaller r' - 
nb color (smaller EW), indicating that they are perhaps 
older galaxies with some dust effects present. 

In the above paragraph, we have outlined different re- 
gions in our color-color plane which would house objects 
from the different scenarios we are examining. We plan 
to use this to determine the likelihood of dust quenching 
of the continuum resulting in a large Lya EW. However, 
the current error bars on individual objects are large 
compared to their distance from the zero-dust line. In 
order to quantify this, we calculated the distance of each 
object from the zero-dust line in units of a. The mean 
distance of the 21 objects we are studying is 1.10 ± 0.48 
a. While it is possible that the dust enhancement of the 
Lya line is present in our sample, due to the large error 
bars we cannot rule out the possibility that all of these 
objects lie on the zero-dust line at a level of ~ 1 a. 

Performing the same exercise with the six stacked 
points with their smaller error bars, the mean distance 
from the zero-dust line is 1.17 ± 0.81 a. However, the 
smaller error bars do allow us to make more of a char- 
acterization of a few of the stacks. We can most likely 
rule out clumpy dust enhancement of the Lya line for 
the low and middle EW stacks, mainly because they lie 
below the zero-dust line, and far from the region in our 
color-color plane where we would expect to find the en- 
hanced line. Also, in the line enhancement scenario, one 
would expect to find extremely high EW Lya lines, and 
these two stacks do not contain the highest 33% EW ob- 



jects in our sample. There still may be dust present in 
these objects, but it would be in a uniform distribution, 
quenching both the continuum and the Lya line. 

5. CONCLUSION 

We have used MMT/Megacam broadband photometry 
in order to study the continuum properties of Lya emit- 
ting galaxies at z^ 4.5. By dividing our objects into 
six different catcrgories and stacking them, we were able 
to derive age and stellar mass estimates for an average 
galaxy in each stack (ignoring dust effects). In all cases, 
the best fit ages were young, with the oldest age com- 
ing from objects with the lowest equivalent width. Even 
those objects had an age of only 200 Myr (40 Myr) for a 
continuous (exponentially decaying) star formation rate. 
As would be expected, the bin with the highest EW had 
the youngest age out of the continuum detected objects 
at 4 Myr (from both SFRs). The youngest age from 
both SFRs came from objects whose continuum flux was 
undetected, mostly because these were detected in the 
narrowband image, and not in the continuum images. 

The derived stellar masses ranged from ~ 2xl0 7 M Q 
for the undetected objects to ~ 2xl0 9 for the lowest 
EW objects. Our derived masses are consistent with 
other mass estimates of similar objects. In the major- 
ity of cases, the masses of our objects were greater than 
those from other studies (Ellis et al. 2001; Gawiser et 
al. 2006; Pirzkal et al. 2007). However, the magni- 
tudes of our individual objects were brighter than those 
of the comparison studies, implying that our objects are 
larger, more massive galaxies than those from the other 
studies. In conclusion, while we have not yet been able 
to definitively determine the likelihood of dust enhance- 
ment of the Lya line, we have been able to shed some 
light on the physical properties of these high-z objects. 
The Lya galaxies in this survey, especially those with 
the largest Lya EW, are some of the youngest and least- 
massive objects in the early universe known to date. 
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